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Abstract 

We have determined lead (Pb) abundances for twelve red giants having [Fe/H] between —2.1 and —1.3, 
and its upper-limits for two lower mctallicity objects, as well as the abundances of lanthanum (La) and 
europium (Eu). The averages of [Pb/Fe] and [Pb/Eu] arc —0.3 and —0.6, respectively, and no clear increase 
of these ratios with increasing metallicity is found. The [La/Eu] values arc only slightly higher than that 
of the r-process component in solar-system material. Combining the previous studies for globular clusters, 
these results suggest little contribution of the s-process to Pb of the field stars studied here, supporting 
the estimate of Pb production by the r-process from the solar-system abundances. 

Key words: Galaxy: halo — nuclear reactions, nucleosynthesis, abundances — stars: abundances — 
stars: Population II 



1. Introduction 

Elements heavier than the iron group are mostly syn- 
thesized by neutron-capture processes, which are distin- 
guished into the rapid (r-) and slow (s-) processes by the 
reaction time scale (e.g. Burbidge et al. 1957). The 
s-process fraction in solar-system material is estimated 
by the classical or stellar s-process models calibrated for 
the nuclei that are yielded only by the s-process (e.g. 
Kappeler et al. 1989, Arlandini et al. 1999). The r- 
process component is derived by subtracting the s-process 
one from the solar-system isotopic abundances. 

The contributions of the s- and r-processes to Galactic 
objects are estimated from the abundance ratios of el- 
ements that represent the two processes (e.g. Ba/Eu, 
La/Eu; Truran et al. 1981; McWilham et al. 1998; Burris 
et al. 2000). Simmerer et al. (2004) investigated the 
La/Eu ratio for a large sample of metal-poor stars, and 
concluded that the contribution of the s-process, with re- 
spect to the r-process, increases with increasing metallic- 
ity, but large scatter of [La/Eu] exists among stars having 
similar metallicity^. 

Another good probe of the s-process at low metallic- 
ity is Pb, one of the heaviest elements including the iso- 
tope with magic numbers for proton and neutron (^°^Pb). 
The Pb abundances have been measured for many carbon- 
enhanced objects, which are believed to be affected by 
mass transfer across binary systems from Asymptotic 
Giant Branch (AGB) stars, and have been providing a 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 
[A/B] = log{NA/NB) - log{NA/NB)Q, and log ea = 
log{NA/Nii) -I- 12 for elements A and B. 



useful constraint on the s-process in low metallicity AGB 
stars. In contrast, measurements of Pb abundances for 
normal (non-carbon-enhanced) stars are quite limited. 
Travaglio et al. (2001) argued enrichment of Pb in the 
Galactic halo, thick and thin disks, reporting the Pb abun- 
dances for several stars. Their sample, however, includes 
a star possibly classified as a CH star, and the results 
for some objects are quite uncertain. Ivans et al. (2006) 
studied the heavy elements, including Pb, for HD 221170 
([Fe/H]= —2.2) in detail. Pb abundances of stars in globu- 
lar clusters were recently determined by Yong et al. (2006) 
and Yong et al. (2007). 

This letter reports our abundance measurements of Pb 
and other elements for bright metal-poor red giants, which 
would belong to the Galactic halo. The results arc com- 
pared with the solar-system s- and r-process abundance 
ratios as well as that of r-process-enhanced, extremely 
metal-poor stars. 

2. Observation, Analysis and Results 

The sample of this study is selected from the list 
of Burris et al. (2000), who determined abundances of 
neutron-capture elements like La and Eu, and that of 
Beers et al. (2000), who studied metallicity and kinemat- 
ics of metal-poor stars. Bright stars {V < 10) were selected 
for our purpose. 

High resolution spectra were obtained for 14 stars given 
in table 1 with the High Dispersion Spectrograph (HDS; 
Noguchi et al. 2002) of the Subaru Telescope during sev- 
eral observing runs including a short program in 2007. 
The resolving power {R = X/6X) of our observations ranges 
between 50,000 and 90,000, and the wavelength coverage 
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Table 1. Object list 



Object 


iV-K)o 


U 


V 


W 


BD+1°2916 


3.025 


85 


-232 


14 


BD+6°648 


2.950 


-165 


-296 


76 


BD+30°2611 


2.979 


-10 


-59 


-284 


HD 3008 


3.201 


64 


-141 


54 


HD 29574(V* HP Eri) 


2.927 


-296 


-213 


-223 


HD 74462 


2.562 


-121 


-284 


70 


HD 141531 


2.912 


-189 


-278 


-67 


HD 204543 


2.716 


-11 


-130 


26 


HD 206739 


2.495 


82 


-107 


-64 


HD 214925 


3.269 


129 


-512 


136 


HD 216143 


2.908 


-329 


-268 


78 


HD 220838 


2.857 


-27 


-76 


22 


HD 221170 


2.433 


-119 


-123 


-22 


HD 235766 


2.396 


-305 


-260 


-4 



0.8 



0.6 



0.4 



0.2 




is 4000-6800 A, 3700-6400 A or 3100-4800 A, depending 
on observing runs. The data reduction was made with the 
standard procedure using the IRAF^ echelle package (e.g. 
Aoki et al. 2005). The signal-to- noise ratios of the spec- 
tra at the wavelength of the Pb I hne used in the analysis 
(4058 A) are 85 (per 0.024 A pixel) or higher. 

Chemical abundance analyses were performed using the 
grid of model atmospheres of Kurucz (1993) with no con- 
vective overshooting (Castelli et al. 1997)^^. The effective 
temperatures (Toff) were estimated from photometric col- 
ors (SIMBAD and 2MASS; Skrutskie et al. 2006) using 
the scale of Alonso, Arribas, Martinez- Roger (1999). The 
reddening was estimated from the dust map of Schlegel, 
Finkbeiner, Davis (1998) and the interstellar Na D line 
absorption (Munari & Zwitter 1997) when available. We 
gave a priority to the {V — K)q given in table 1 in the 
effective temperature determination, while other colors 
were also investigated for comparison purposes. The sur- 
face gravity (g), micro-turbulence (titurb), and metallicity 
([Fe/H]) are determined by the LTE analysis for Fe I and 
Fe II lines, as to obtain no correlation of Fe abundances 
with line strengths and ionization stages. For the analysis 
of most objects, we adopted the line list of Ivans et al. 
(2006). We added Fe I and II lines from O 'Brian et al. 
(1991), Bard, Kock, and Kock (1991), and Fuhr, Martin, 
Wiese (1988) for relatively metal-rich stars, for which only 
spectra of the UV-blue range are available. Although a 
slightly lower logg value than zero is derived for HD 3008 
and HD 214925, we adopted \ogg = 0.0 for these two stars. 
This little affects the derived abundance ratios of neutron- 
capture elements. The stellar parameters adopted for the 
following analyses are given in table 2. The details of 
the determinations of stellar parameters will be reported 
in our future paper (Aoki et al., in preparation). The 



IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation, 
http: / / wwwuser . oat . ts . astro, it / castelli / grids . html 
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Fig. 1. Comparisons of synthetic spectra for the region in- 
cluding the Pb I 4058 A line with the observed spectrum. 
The Pb abundances assumed in the calculations are loge (Pb) 
= +0.43 ±0.10 (BD+30°2611), 0.00 ±0.15 (HD 204543) and 
-0.1,-0.3 (HD 216143) for solid lines, and -oo for dotted 
lines. The contamination of the CH line at 4057.7 A is not 
severe in the spectra of our non-carbon-enhanced stars. 



carbon abundances are estimated from the CH band at 
4223 A. The [C/Fe] of our sample ranges between —0.2 
and —0.9, indicating these stars are not carbon-enhanced 
objects (table 3). 

The Pb abundances were determined by fitting syn- 
thetic spectra to the observed Pb I 4058 A line. The 
line list of Simons et al. (1989) including hyperfine split- 
ting and isotope shifts, which was used by Aoki et al. 
(2002), was adopted, assuming the solar Pb isotope ra- 
tios. We note that the derived abundance is not sensitive 
to this assumption. Examples of the observed and syn- 
thetic spectra are shown in figure 1. 

The La and Eu abundances are determined using the 
line lists of Lawler, Bonvallet, Sneden (2001a) and Lawler 
et al. (2001b), respectively (see also the appendix of Ivans 
et al. 2006), including hyperfine splitting. The Eu iso- 
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tope ratio (^^^Eu/^^'^Eu) is assumed to be unity. The 
La abundance is determined from more than ten lines for 
all objects, while the Eu lines used are dependent on the 
targets, because of the differences of the wavelength cov- 
erage and the line strengths. The Eu abundances are de- 
termined from lines in the red region (6049 A, 6436 A, 
and/or 6645 A) for most stars. Two lines in the blue re- 
gion (4129 A and 4205 A) are used for objects for which 
only blue spectra are available. We confirmed that the 
Eu abundances derived from the blue and red lines agree 
within 0.13 dex for HD 220838 and HD 221170. The re- 
sults of the La, Eu and Pb abundances (loge values) are 
given in table 2, and the abundance ratios are given in 
table 3. 

The random errors of the derived abundances are es- 
timated as aN~^l'^ , where a is the standard deviation of 
the derived abundances from individual lines and N is the 
number of lines used. The cr values for Fe and La lines are 
~ 0.1 dex or smaller. We adopted cr of Fe I lines as that of 
Eu, for which five lines or less are available for abundance 
analyses. The random errors of the Pb abundances are 
estimated from the fitting of the Pb I line, including the 
uncertainty of continuum placement. For stars showing a 
clear Pb I feature (the top panel of figure 1), 0.10 dex is 
adopted as the fitting error, while 0.15 dex is adopted for 
stars having weaker features (the middle panel of figure 1). 
No clear feature is found for HD 216143 (bottom panel of 
figure 1) and BD-|-6°648, for which we determined upper 
limits of Pb abundances. 

Errors due to the uncertainty of atmospheric pa- 
rameters are estimated by calculating abundances 
changing atmospheric parameters by A(rcff) = 100 K, 
A(log5)=0.3 dex, A([Fe/H])=0.2 dex, and A(vturb) = 
0.2 km s-i for BD-f30°2611. These errors arc added, 
in quadrature, to the random errors estimated above, and 
are given in table 2. The effects of the changes of at- 
mospheric parameters on the abundance ratios of Pb/La, 
Pb/Eu, and La/Eu are also calculated, and are included 
in the errors given in table 3. We note that the Pb abun- 
dance derived from the neutral species is lower if higher 
gravity is adopted, while the sign is opposite for La and 
Eu abundances from ionized species. This results in larger 
errors in [Pb/La] and [Pb/Eu] than those in [La/Eu]. 

A comparison of our abundance results for HD 221170 
with those of Ivans et al. (2006) shows 0.1-0.2 dex differ- 
ences in the abundances of Fe, La, Eu and Pb, which 
are, at least qualitatively, explained by the differences 
of adopted atmospheric parameters (ATes =90 K and 
Alog5 =0.5 dex). The abundances of HD 141531 are 
determined by the present work for the same spectrum 
as used by Yong et al. (2006). The results agree within 
0.05 dex between the two measurements, adopting almost 
the same stellar parameters. 

3. Discussion and concluding remarks 

The top panel of figure 2 shows [Pb/Fe] as a function 
of [Fe/H] for the field stars of our sample and for glob- 
ular cluster stars. The Pb of our sample is underabun- 
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Fig. 2. Abundance ratios of tlie neutron capture elements 
La, Eu, and Pb as functions of [Fe/H], Circles indicate our 
results. Open circles mean the result for the two objects for 
which only an upper limit of the Pb abundance is determined. 
Open triangles are the results for the globular clusters M4, 
M5, M13, and NGC 6752 by Yong et al. (2006) and Yong et 
al. (2007). The dotted lines in the lower three panels show the 
abundance ratios of the solar-system r-process components 
(Simmerer et al. 2004). 

dant, and the abundance scatter is small: the average of 
[Pb/Fe] is —0.27 and the standard deviation is 0.12 dex 
(table 4). This result indicates no rapid increase of Pb 
in this metallicity range in our sample. The open trian- 
gles in this panel present the averages of Pb abundance 
ratios of stars in the globular clusters M 4, M 5, M 13 and 
NGC 6752 (Yong et al. 2006; Yong et al. 2007). While 
the [Pb/Fe] of three clusters agrees with that found for 
our field stars, the value of M 4 ([Fe/H]= -1.23) is higher 
by 0.5 dex. The high Pb abundance ratio should be the 
result of the large s-process contribution to this cluster 
(Ivans et al. 1999). 
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In order to investigate the origins of Pb in the field 
stars, we depict the abundance ratios of Pb/La, Pb/Eu 
and La/Eu in the lower three panels of figure 2. The 
dotted lines in these panels mean the abundance ratios of 
the r-process component in solar-system material (table 4) 
determined by Simmerer et al. (2004). The Pb/Eu and 
La/Eu of the s-process component are much higher than 
those of the r-process component. Hence, the agreement 
of [Pb/Eu] of the field stars with that of the solar-system 
r-process component implies that the r-process is the dom- 
inant source of Pb in our sample, and the contribution of 
the s-process is small if any. 

The little contribution of the s-process to our sample 
is also suggested by the La/Eu ratios. The La and Eu 
abundances of a large sample of metal-poor stars were in- 
vestigated by Simmerer et al. (2004). Seven stars among 
our objects are included in their sample, and the [La/Eu] 
values measured by the two studies agree well. However, 
their sample also includes stars having relatively high 
La/Eu ratios ([La/Fc]~ —0.2) in this metallicity range, 
and they concluded the existence of a spread in the level 
of s-process enrichment. Measurements of Pb for the stars 
having higher La/Eu ratios are desirable. It should be 
noted, however, that most of such stars in their sample 
have relatively high temperature, and measurements of 
Pb abundances are relatively difficult. 

Simmerer et al. (2004) argued a correlation between the 
La/Eu ratio and the kinematic properties of the objects, 
suggesting that "high velocity stars" having high U and 
W, or low V, show low La/Eu ratios. The kinematics 
data of our objects taken from Beers et al. (2000) are 
given in table 1. They show the characteristics of "high 
velocity stars", with a possible exception of HD 220838, 
supporting the arguments of Simmerer ct al. (2004). 

The above discussion is based on the r-process com- 
ponent in solar-system material estimated by Simmerer 
et al. (2004). However, the Pb productions by the s- 
and r-processes are still very uncertain. Indeed, the 
detailed abundance studies for r-process-enhanced, ex- 
tremely metal-poor stars reported significantly lower Pb 
abundance ratios (Plez et al. 2004; Frebel et al. 2007). 
For instance, CS 31082-001 has [Pb/Eu] of -1.28 (Plez 
et al. 2004; Hill ct al. 2002), 0.63 dex lower than that 
of the solar-system r-process component (table 4). If 
we adopt this as the abundance ratio of the r-process 
contribution (assuming Eu is entirely the r-process ori- 
gin), only 20% of Pb of our sample (<[Pb/Eu]>= -0.59) 
is from the r-process and the remaining 80% is the s- 
process origin. In contrast, the average of [La/Eu] of 
our sample (<[La/Eu]>=: —0.41) is only slightly higher 
than that of CS 31082-001 ([La/Eu]= -0.52). If simi- 
lar estimate to Pb is applied, 22% of La originates from 
the s-process. These estimates results in the Pb/La ra- 
tios of the "s-process component" of the field stars to be 
[Pb/La]s =<[Pb/La]> +log(0.80/0.22) = 0.38. 

Thus, an alternative interpretation for our result is that 
the r-process yield of Pb is as low as that found in the r- 
process-enhanccd extremely metal-poor stars, while the 
Pb production with respect to La by the s-process at low 



metallicity is quite efficient. An efficient synthesis of ^"'^Pb 
at low metallicity is predicted by the s-process calcula- 
tions for AGB models with ^^C pocket (e.g. Busso et al. 
1999). However, the decrease of Pb/La ratios (possibly 
by one dex) with increasing metallicity, which is expected 
from such models, is not found in our sample. Another 
difficulty of this interpretation is the abundance ratios of 
M 4 stars. While the high La and Pb abundances of this 
cluster indicate large contributions of the s-process, the 
< [Pb/La] > of this cluster (table 4) is significantly lower 
than the [Pb/La]s estimated above. That is, assuming 
large contributions of the s-process to our field stars re- 
quests an unnatural assumption that the s-process yields 
are significantly different between field stars and globular 
cluster objects. 

Our conclusion here is that Pb, as well as other neutron- 
capture elements, in a majority of halo stars in the metal- 
licity range of [Fe/H]< —1.3 primarily originated from 
the r-process. This supports the Pb production by the 
r-process estimated from solar-system abundances (e.g. 
Burris et al. 2000; Simmerer 2004), and suggests that the 
Pb abundance of r-process-enhanced, extremely metal- 
poor stars (e.g. CS 31082-001) is anomalously low. Such 
diversity of Pb production by the r-process is recently ar- 
gued by Wanajo (2007). However, given the fact that 
some stars having high La/Eu ratios are known in this 
metallicity range, further measurements of Pb abundances 
for a larger sample are strongly desired. 
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Table 3. Elemental Abundances Ratios 
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.93 


— 


.51 


—0 42 


0. 


22 


nil 


0. 


26 


— u.oo 


0. 


.27 


— 


.42 


0. 


.21 


HD 29574 


-1, 


.81 


-0, 


.83 


-0.34 


0. 


19 


—0.30 


0. 


,22 


—0.73 


0. 


.25 


-0, 


.43 


0. 


.20 


HD 74462 


-1, 


.34 


-0, 


.40 


-0.31 


0. 


,26 


-0.25 


0. 


,29 


-0.59 


0. 


.31 


-0, 


.34 


0, 


.20 


HD 141531 


-1, 


.61 


-0. 


.58 


-0.34 


0. 


19 


-0.18 


0. 


,22 


-0.56 


0. 


.25 


-0, 


.38 


0. 


.20 


HD 204543 


-1, 


.81 


-0. 


.63 


-0.19 


0. 


,22 


0.00 


0. 


,25 


-0.18 


0. 


.27 


-0, 


.18 


0, 


.20 


HD 206739 


-1, 


.59 


"0. 


.25 


-0.21 


0. 


,26 


-0.17 


0. 


,29 


-0.56 


0. 


.31 


-0, 


.39 


0, 


.20 


HD 214925 


-2, 


.03 


-0, 


.81 


-0.47 


0. 


,22 


-0.42 


0. 


,25 


-0.89 


0. 


.27 


-0, 


.47 


0, 


.20 


HD 216143 


-2, 


.27 


-0. 


.37 


< -0.69 






< 0.33 






< -0.34 






-0, 


.67 


0, 


.20 


HD 220838 


-1, 


.75 


-0. 


.49 


-0.20 


0. 


19 


0.03 


0. 


,22 


-0.50 


0. 


.25 


-0, 


.53 


0, 


.20 


HD 221170 


-1, 


.98 


-0. 


.61 


-0.22 


0. 


,26 


-0.44 


0. 


,29 


-0.95 


0. 


.31 


-0, 


.51 


0, 


.20 


HD 235766 


-1, 


.88 


-0, 


.66 


-0.02 


0. 


,26 


-0.08 


0. 


,29 


-0.52 


0. 


.31 


-0, 


.44 


0, 


.20 



Table 4. Abundance ratios of neutron-capture elements 



Object 


Ref.* 


[Pb/Fe] 


[Pb/La] 


[Pb/Eu] 


[La/Eu] 


field stars (12 objects)^' 


1 


-0.27 (0.12) 


-0.18 (0.15) 


-0.59 (0.20) 


-0.41 (0.09) 


M 4 stars (12 objects)!' 


2 


0.30 (0.07) 


-0.18 (0.06) 


-0.09 (0.09) 


0.08 (0.06) 


s-process (solar-system) 


3 




0.022 


1.471 


1.449 


r-process (solar-system) 


3 




-0.057 


-0.650 


-0.593 


CS 31082-001 


4 


0.35 


-0.76 


-1.28 


-0.52 



* References: (l)This work; (2)Yong et al. (2007); (3)Simmerer et al. (2004); (4)Hill et al. (2002) and (Plez et al. 2004). 



t Averages and standard deviations of abundance ratios. 



